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^ ' Abstract 

^ . Higher order corrections to the Bahtsky-Kovchegov equation have been 

^ I estimated by introducing a rapidity veto which forbids subsequent emissions to 

be very close in rapidity and is known to mimic higher order corrections to the 
Unear BFKL equation. The rapidity veto constraint has been first introduced 
using analytical arguments obtaining a power growth with energy, (5s(Y) ~ e^^, 
of the saturation scale of A ~ 0.45. Then a numerical analysis for the non- 
linear Balitsky-Kovchegov equation has been carried out for phenomenological 
rapidities: when a veto of about two units of rapidity is introduced for a fixed 
value of the coupling constant of = 0.2 the saturation scale A decreases from 
~ 0.6 to ~ 0.3, and when running coupling effects are taken into account it 
decreases from ~ 0.4 to ~ 0.3. 
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1 Introduction 

The high energy behavior of a parton system can be associated to the Bahtsky-Fadin- 
Kuraev-Lipatov (BFKL) dynamics pp. At leading order (LO) the BFKL equation 
resums contributions of the form (a^Y)", with Y ~ Ins being a rapidity variable. 
Linear evolution gives rise to a Pomeron-like behavior of the scattering amplitudes 
with an intercept u^^^^ ~ 0.5. This power growth of the amplitude with energy 
violates s-channel unitarity at rapidities of the order of Y ~ l/cts In l/a^ 

A theoretical possibility for the high energy growth of the amplitudes to be mod- 
ified in a way consistent with unitarity is the idea of parton density saturation 
which accounts for the possibility of parton annihilation, an essentially nonlinear ef- 
fect. Present theoretical understanding views a system of saturated partons as a new 
state of matter called Color Glass Condensate (CGC) (see e.g. Ref. ^ and references 
therein). Saturation signals have been intensively searched experimentally mostly 
at HERA and more recently at RHIC. Although there are hints of parton density 
saturation in the present data, a clean and ultimate signal has not been found so far. 

The fundamental quantity characterizing the transition to the saturation regime 
is the so-called "saturation scale", Qs{Y)- The determination of the rate of growth 
with rapidity of this saturation scale could be of a large importance for, e.g., structure 
function extrapolations from HERA to LHC kinematics. In the context of the well- 
known saturation model of Golec-Biernat and Wusthoff (GBW) [S], the saturation 
scale grows exponentially as Qs(Y) ~ exp(A Y/2) with A ~ 0.29. 

Within the LO approximation a theoretical tool with solid grounds in perturba- 
tive QCD suitable to study saturation phenomena at high energies is the Balitsky- 
Kovchegov nonlinear evolution equation (BKe) ^. This equation reads 



with as = otsNc/ vr. In the color dipole approach to high energy scattering the function 
A^(rx,Y,6) stands for the imaginary part of the amplitude for a dipole of size 
elastically scattered at an impact parameter h. In this paper the impact parameter 
dependence of the amplitude will be neglected, considering, in this way, a target of 
infinite size, p is an ultraviolet cutoff needed to regularise the integral which does not 
appear in physical quantities. 

The physical content of Eq. is that of a dipole of size xio which decays into 
two other dipoles of size Xi2 and xo2 with a decay probability given by the wave 

^Eq. was originally proposed by Gribov, Levin and Ryskin ^ in momentum space and derived 
in the double logarithmic approximation by Mueller and Qiu [3| . In the leading In 1 /x approximation 
it was obtained by Balitsky using a Wilson Loop Operator Expansion 7 . In the form presented 
in Eq. ^ it was obtained by Kovchegov 8; using the color dipole approach 9 to high energy 
scattering in QCD. This equation was also obtained by summation of BFKL Pomeron fan diagrams 
by Braun TU^ and, more recently, Bartels, Lipatov and Vacca In the framework of the Color 
Glass Condensate it was obtained by lancu, Leonidov and McLerran ^21 




L- [2Ar(xo2, Y) - Ar(xoi, Y) - Ar(xo2, Y)iV(xi2, Y)] , (1) 
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9 X m 

function |\&| = 2 ■ These two dipoles can then interact with the target and, 
^02 ^12 

with a certain probabihty, they can do so simultaneously, a possibility accounted for 
by the non-linear term in Eq. (^J. The linear part in Eq. (PJ corresponds to the LO 
BFKL equation describing the evolution with energy of the multiplicity of the fixed 
size color dipoles. The BKe has been studied both analytically jTSl El El Ej and 
numerically [inilIIllIHllinil2niinil2ai23l21|- Phenomenologically the BKe provides 
a good description of DIS data from HERA [13 123 ESI l2Z|- It is worth pointing 
out that the linear part of the BKe is obtained in the leading soft gluon emission 
approximation keeping the strong coupling fixed and that the large Nc limit is used 
in order to write the nonlinear term as a product of two functions N. This limit 
is at the basis of the color dipole picture and, to a large extent, it corresponds to 
a mean field theory without dipole correlations. The equation also neglects target 
correlations, an assumption which might be valid for asymptotically heavy nuclei but 
not for protons or a realistic nucleus target. 

It would be very interesting to go beyond the original BKe and relax some of the 
underlying assumptions outlined above. At present there is a large activity in this 
direction. Regarding the contribution of the Nc corrections they can be estimated to 
be up to 15% ^21,. In this publication we would like to focus on the higher order 
corrections which are relevant, in particular, for phenomenological applications. 

In principle, unitarity corrections based on LO estimates are expected to be im- 
portant at rapidities of the order Y ~ l/a^ In l/a^, parametrically earlier than the 
next-to-leading (NLO) corrections which set in at Y ~ l/tt^- It is also known that 
the NLO corrections to the linear BFKL equation significantly decrease the Pomeron 
intercept thus postponing the arrival of unitarity corrections to higher rapidities. 

A complete nonlinear equation at NLO has not been derived yet. In the conven- 
tional approach based on s-channel unitarity, the forward BFKL kernel is known at 
NLO [23 I2n|- A nonlinear evolution needs the knowledge of the non-forward ker- 
nel |HD| together with the NLO impact factor [Hll IS2l IHHj, both calculations being 
currently under investigation. However, a NLO study of the triple Pomeron vertex 
entering the BKe has not been initiated yet. So far, the only exact result which has 
been reported is due to Balitsky and Belitsky [34j who have been able to compute a 
single NLO contribution with maximal nonlinearity, the A^^ term. 

There have been some attempts to get insight about saturation at NLO using 
approximate methods. Triantafyllopoulos jHSj has considered the renormalisation 
group improved NLO BFKL equation with the presence of a saturation boundary. 
His results show a decrease in the saturation scale growth as a function of rapidity 
towards the value A ~ 0.3 observed experimentally. A similar type of study based 
on the NLO BFKL has been recently reported in !HB]j. 

In this work we propose a new approach for the study of saturation effects includ- 
ing NLO corrections. We will introduce a constraint in the rapidity of the emitted 
gluons in the BKe, a so-called "rapidity veto" [SZl IHH] which, for the linear part of 
the equation, is known to reproduce the bulk of the NLO corrections. In the next 
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section we revise how to introduce a veto in rapidity in the hnear BFKL equation 
and apply this constraint to obtain an estimate of the saturation scale as a function 
of the veto. In Section El we apply the method of rapidity veto to the BKe and study 
its influence on the energy growth of the saturation scale. In the last Section of this 
work we present our summary. 



2 The rapidity veto in BFKL and the saturation 
line 

In the following the introduction of a rapidity veto as in Ref. EH| will be shown. 
To impose the constraint that subsequent gluon emissions are separated by some 
minimum interval in rapidity, rj, can be done writing the LO BFKL equation as an 
integral equation in rapidity, i.e. 

/(Y,7) = JdY'e{Y-Y'-r^)asX{l)f{Y'n), (2) 

where 7 corresponds to a Mellin transform in transverse momentum space and xil) = 
2\1/(1) — \l/(7) — — 7) is the eigenvalue of the LO kernel. To go to the representation 
in the u plane we use the transformation 

fM = j dYe-^^ f{Y,^). (3) 

We can now introduce this in the original equation to obtain 

Ul) = «.X(7) /ciY7(Y',7) / ciYe-^Y, (4) 

J JY'+r) 

therefore the effect of imposing the veto on the LO BFKL equation leads to an 
eigenvalue which is determined by the solution to 

UJ = a.x(7)e-''^ (5) 

It is worth noting that the solution to this equation respects the structure of a max- 
imum at 7 = ^ + i 1/ for 1/ ~ so that this region dominates at high energies. This is 
highlighted in Fig. ^ where the maxima are shown revealing how the original value of 
the Pomeron intercept decreases from about 0.5 to about 0.3 for a value of the veto 
of two units of rapidity. This is in agreement with other predictions from studies of 
the NLO gluon Green's function [HH] . 

Now we calculate the influence of this veto in rapidity on the saturation scale. In 
the case of forward scattering the amplitude for a dipole of size 1/Q on a dipole of 
size l/Qo can be written as 

N{Q,Qo,Y) = j d^ j dujNo{l)e^v{-lL + uY) (6) 

J J uj-asX (7) 
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Figure 1: Dependence of the LO eigenvalue of the kernel on the veto upon u 



with L = InQ^/Qo- '^^^ ^^to in rapidity is easily introduced via a modified kernel as 
in Eq. ®: 



Ar(Q,Qo,Y,r7) 



(i7 / duAfo (7) exp (— 7-L + luY) 



Lu - asX (7) e 



(7) 



The saturation line, Ls = lnQ'j.{Y)/Ql with (5s(Y = 0) = Ql, can be defined as that 
with a stationary exponent: 

-7i^.+^(7,^)Y = 0, (8) 

where the introduction of the veto enforces 

w(7,r/) = a,x(7)exp(-r/^(7,r/)). (9) 

At high energies the dominant region is that in the intersection with the saddle point 
7 



Y 



0. 



(10) 



7=7 



The solution to this system of equations provides an implicit equation for 7: 
^7-1 - «..x(7)exp 1-^7 . 



(11) 



0.71 




0.62 I 1 1 1 1 1 1 

0.5 1 1.5 2 2.5 3 

n 

Figure 2: The solution to Eq. (fTTjl as a function of the veto in rapidity, rj. 

Consequently, when the rapidity veto is imposed it develops a dependence on the agT] 
product, 7 = 7 (Os?/), and the saturation line reads now 

Ls = as^Yexp(^l-^^^^X{as,v)Y. (12) 

For a value of = 0.2 in Fig. |21we have plotted the dependence of 7 on the rapidity 
veto rj. The effect of the constraint in rapidity is to increase the value of this critical 7 
from the well known 7 ~ 0.63 (a value which is obtained for rj —>■ and/or in the limit 
of as 0) to about 0.69 for rj ~ 2.5. This is in agreement with the recent results for 
this quantity of Ref. [SE] where a resummed NLL BFKL equation was under study. 

The evolution in energy is determined by the dLg/dV = A(q;s,?7) derivative. This 
A term is calculated in Fig. El where it can be seen how the effect of the rapidity 
constraint is to delay the onset of the saturation line by means of decreasing A at 
larger values of the veto. For zero veto it corresponds to the usual value of this linear 
coefficient of A ~ 0.93, which can be calculated from Eq. (fTT|l setting t] = 0, a.^ = 0.2 
and then introducing the result for 7 in Eq. (fT^ reading A = asX(7)/7- At a rapidity 
constraint of rj = 2.5 we obtain A ~ 0.45, in agreement with Ref. [3^ and larger than 
that calculated in Ref. [53] . 

These results have been produced in the case of linear evolution imposing a con- 
straint which reproduces higher order corrections and estimating the position of the 
saturation line. In the following section we introduce the veto in rapidities directly 
in the non-linear BKe and study the consequences of this constraint on the evolution 
for phenomenological rapidities. We will see that the power growth of the saturation 
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Figure 3: The saturation power (Eq. ((121)) a function of the veto in rapidity, rj. 

scale is slower, even for zero veto, mainly due to preasymptotic effects. We will then 
show how the effect of the veto in rapidity is very similar to that found in this section: 
the effective power decreases as the veto is larger. 

3 The rapidity veto in the BKe 
3.1 The fixed coupling case 

To proceed with the numerical analysis it is convenient to write the BKe for fixed 
coupling in the integral form 

iV(xoi,Y) = iV(xoi,Yo) + as f dY' [ x 

JYo J p 2,71 Xo2 Xf 2 

[2iV(xo2,Y') - iV(xoi,Y') - iV(xo2,Y')iV(xi2,Y')], (13) 

with the initial condition being defined at the rapidity Yq. As it was said above 
the rapidity veto prevents two emissions from being emitted close to each other in 
rapidity space. In this work we impose this veto in both the linear and non-linear 
parts of the BKe. In this way we maintain the locality of the recombination process, 
corresponding to the quadratic term. Hence, the new equation simply reads 

iV(xoi,Y) = iV(xoi,Yo) + as f ~' dY' f J"^ x 

[2iV(xo2,Y') - iV(xoi,Y') - iV(xo2,Y')iV(xi2,Y')]. (14) 
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In the numerical implementation we work with the differential form of this equa- 
tion with veto, i.e., 



which highlights the non-locality in rapidity after the constraint has been imposed. 
Expanding Eq. (|T3j) in t] (we assume 77 -C Y) one can easily verify that the corrections 
introduced via veto are proportional to al and thus of the NLO. 

In the context of the BKe the imposition of the veto has the consequence of a 
small fluctuation above = 1 when the function approaches the unitarity bound. 
These small fluctuations do not grow with rapidity. Technically this small violation 
of unitarity has its origin in the fact that the evolution is not stopped at A^(Y) = 1 
since the derivative is computed at a retarded rapidity Y — 77, see Eq. fll5|l . This small 
instability is not surprising since our approach is only an estimate of the NLO correc- 
tions. To fully preserve unitarity possibly we would have to introduce a correlation 
between coordinates and rapidity. Connecting with this point it is worth noticing 
that a generalization of the BKe proposed in Ref. jlO] and amounting to having extra 
1 — A^(xoi, Y) factor in front of the evolution kernel would respect unitarity even in 
the presence of a rapidity veto. In order to study the behavior of the saturation scale, 
our analysis will be centered around the transition region ~ 0.5, which is a region 
not affected by the above issue. 

Eq. ()15|) will be solved numerically for x < Xq = 0.01, which corresponds to 
rapidities above Yq ~ 4.65. With this goal in mind we need to specify the initial 
conditions to the non-linear equation which, in principle, should be fitted to experi- 
mental data. Motivated by the phenomenological accuracy of the results in Ref. jl7j . 
the same initial conditions as those in that reference are used in the present work. 
These conditions were fitted to low x F2 data for the BKe with no veto and read 



Here is taken to be LO running at the scale and xG'^'^^^ is the LO CTEQ6 
gluon distribution also computed at the scale The initial condition in (fTBj) is 

smoothly extrapolated to = 1 at very large distances using the method proposed 
in Ref. j2E] and implemented in Ref. The parameter R stands for the effective 
proton size, = 3.1GeV~^, an output of the fit performed in Ref. [Ej. For the 
numerical implementation of the veto it is necessary to generate the initial conditions 
in a band of width t] between Yq and Yq — rj. The reason for this becomes clear if 
we think the evolution in rapidity Y as a process in which N{r, Y) is the input for 
the next step which will give us N{r,Y + 6Y). In the non zero veto case in order to 
compute A^(r, Y + 6Y) one would need to consider N{r,Y — rj + 6Y). Hence, since 
we start at Yq we need a band of initial conditions that will span the space between 
Yq — ?7 and Yq. As there is no evolution in energy along that band we assume the 
initial conditions to be independent of Y and equal to N{r, Yq) on the band. 




[2Ar(xo2,Y-r/) - Ar(xoi,Y-r/) - Ar(xo2, Y - r^) Ar(xi2, Y - 77)] , (15) 



N{r,Yo) = 1 - exp(-a,C^r2xG^T^Q/(7ri?2)). 



(16) 
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In Fig. m the first result for tlie solution to the BKe is shown. There it can 
be seen how the amplitude starts at zero for small transverse distances r (color 
transparency) to reach the saturation regime ~ 1 at r about 4 GeV~^. This trend 
is general for all values of the rapidity veto but, as a new feature, we observe how 
the arrival of the saturation of the amplitude is delayed as the veto increases. This 
first plot was done for a rapidity of 10 and a fixed coupling of = 0.2. In Fig. El 




Figure 4: Solution to BKe for different values of the veto as a function of r, for Y = 10 
and as = 0.2. 

we highlight how in a region of larger rapidity, Y = 14 saturation arises earlier in r 
for the set of initial conditions we have chosen. In this case of larger center-of-mass 
energy the effect of the veto is more dramatic, considerably delaying the onset of 
saturation. 

The major effect of the rapidity constraint comes when studying how the BK 
amplitude evolves with energy. As the intercept in the linear part is significantly 
reduced when higher order corrections are taken into account the saturation of the 
amplitude comes also later in rapidities. To illustrate this point we plot Fig. 
where we have chosen a typical value of r = 0.75 GeV~^. Once again our estimated 
NLO corrections do delay the onset of saturation. To make this statement more 
quantitative we now proceed to study the saturation scale and its dependence with 
energy. Different definitions of the saturation scale can be associated with the step 
like function A^ [211 ■ They might lead to different normalizations although the energy 
dependence is qualitatively unique. For the sake of simphcity here we adopt the choice 
proposed in Ref. 0^ where it was taken at a point where A^ reaches half: 
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Figure 5: Solution to BKe for different values of the veto as a function of r, for Y = 14 
and as — 0.2 
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(17) 



with Vs = 2/Qs. For phenomenological applications the behavior of InQ^ with rapid- 
ity Y can be fitted hnearly as in Eq. (fT2|) . i.e., 

Ql = Qle'\ (18) 

The numerical analysis of the rapidity dependence of this saturation scale is carried 
out in Fig. [7| This plot refiects very clearly how saturation tends to appear later in 
rapidity, in particular, for the veto which reproduces the NLO intercept, r/ ~ 2. We 
have also performed a linear fit to estimate the linear power of Eq. p2|) . A, this fit was 
done for phenomenological rapidities between Y = 10 and Y = 15 so the expected 
value of the growth cannot be as large as in Fig. IHldue to preasymptotic effects ^. In 
fact, A is much smaller at the beginning of the evolution for rapidities up to Y ~ 10. 
In Fig. ISlwe observe the transition power being of the order of 0.65 already at zero 




Figure 7: Rapidity dependence of the saturation scale for fixed and different vetoes. 

rj. This value is smaller than the equivalent obtained in Sec. |21 Remarkably, the 
dependence on the rapidity veto is of the same functional form as in Fig. IHl with A 
reaching ~ 0.31 at a veto of ?7 = 2.5. 

The main conclusion of this section is to confirm the delay in energy space of the 
arrival of saturation when estimated higher order corrections are introduced in the 

^The preasymptotic effects in saturation scale have been studied in Refs. |41l 1151 llHj . The 
numerical size of these terms is large at lower rapidities with the asymptotic values only reached at 
very large rapidities of the order of Y ~ 100 (see Ref. 1^ for a similar discussion). 
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Figure 8: t] dependence of A with running and fixed = 0.2. 



BKe with fixed couphng constant. The numerical results when we introduce the veto 
on the BKe are consistent with those obtained from a more analytical approach in 
Sec. El based on the BFKL equation with a constraint in rapidity. In the following 
section the effect of this veto will be studied for phenomenological rapidities and 
running the QCD coupling. 



3.2 The running coupling case 

The BKe was originally derived for constant a^, the introduction of the running is 
part of the NLO corrections. At present the use of the running in the BKe can 
only be done by modeling. In previous analysis introducing running seems to be 
phenomenologically favored by the data ^^I- This is because the effect of running a.^ 
is to bring A down to about 0.3 in the phenomenologically relevant region of x > 10^^. 
The main concern of this section is to study the stability of this value of A ~ 0.3 when 
the veto is imposed on top of the running. 

Similar to the initial conditions, as is taken at the leading order running with 
respect to the external scale At large distances as is frozen at the value ~ 0.5. 

We have checked that our results are not sensitive to variations of this value. 

Let us start with Fig.lHlwhere we again show the region of small r for the amplitude 
as calculated from the BKe introducing the rapidity constraint, this time running the 
coupling. As previously found the effect of the higher order corrections is to delay 
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saturation. The rapidity chosen for this plot is 14. It is worth pointing out that the 
effect of the veto is reduced if we compare Fig. IHl to Fig. we will go back to this 
point soon below. 




What about the energy dependence of the saturation scale? The answer to this 
question is plotted in Fig. El where the saturation scale is shown as a function of 
rapidity. The usual delay of the onset of saturation can be again observed although 
the effect of the estimated higher order corrections is smaller than for the case of 
fixed coupling in Fig. [7| To make this more explicit we calculate the dependence of 
the A power in Eq. (|T2|) as a function of the rapidity veto t] in the case of running 
coupling ^ (Fig. (jH}). The main conclusion is that the effect of the veto is not so big in 
the running coupling case, taking A from ~ 0.37 for r] = to about 0.27 for t] = 2.5. 
This variation is much smaller than for the fixed coupling case. 

The calculations in this section teach us that once the coupling is allowed to 
run the influence of other higher order corrections is diminished. The prediction for 
the growth of the saturation scale with energy remains of the order of A ~ 0.3 for 
phenomenological energies independently of the rapidity veto. 

As a final remark, we have checked that the solutions to the BKe for both fixed and 

■^It is known from analytic studies that, contrary to the fixed coupling case, for running coupling 
InQs VY- However, again due to large preasymptotic corrections 24 , InQ^ can be fitted linearly 
in rapidity for a limited range in Y relevant for phenomenology. 
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Figure 10: Y dependence of the saturation scale for running a^- 

running cases exhibit the geometrical scaling property |12j. Namely, the amplitude 
N{r,Y) is a function of the product r = rQs(Y), that is N{t). Scaling holds in the 
saturation domain r > 2 and extends to a much broader region r << 2. For the BKe 
with rj = scaling was shown in Refs. [THl 1211 dl • We have found that the scaling is 
still present after introducing a non-zero veto, suggesting the full NLO amplitude to 
be also a scaling function. 

4 Conclusions 

In this work higher order corrections to the Balitsky-Kovchegov equation have been 
estimated. This estimate has been based on the introduction of a so-called "rapidity 
veto", which forbids two emissions to be very close in rapidity. It is known from 
Ref. P?! IHH] that the introduction of a veto, r], of ~ 2.5 units of rapidity mimics 
the Pomeron intercept predicted by other resummations of the NLO BFKL equation. 
We have estimated these higher order corrections first using analytical arguments 
imposing the rapidity veto on the LO BFKL equation obtaining a power growth of 
the saturation scale of A ^ 0.45, for = 0.2, consistent with that calculated in 
Ref. and larger than that obtained in Ref. [33], our main result being Eq. pi|) . 
Eq. ()12|) and shown in Fig. El 

We have then pursued a numerical analysis of the introduction of the rapidity 
constraint in the full Balitsky-Kovchegov equation for phenomenological rapidities, 
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without using asymptotic arguments. For a fixed coupling constant of 0.2 we observe 
tliat tlie power A decreases from ~ 0.65 for zero veto to ~ 0.31 for rj = 2.5, with a 
dependence on the veto very similar to the previous analytical study, the main result 
plotted in Fig. |H1 When running coupling effects are also taken into account the effect 
of imposing the veto is not so important taking A at no veto from ~ 0.37 to ~ 0.27 
for r] = 2.5 (Fig. El). 

As follows from the present analysis the running coupling effects account for the 
bulk of the NLO corrections to the BKe (see Fig. Given that, as shown in |24j . 
the dipole amplitude is not sensitive to the way the running is introduced, this sug- 
gests that phenomeno logical analysis including running coupling effects, as in fTf], do 
provide reliable predictions. 
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